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The non-dipolar magnetic fields of accreting T Tauri stars 



S. G. Gregory 



oo 
O 
O 
(N 

5 1 

6 

& 



> 

00 

in 

q 

o 
oo 
o 



% 



l * S. P. Matt 2 , J.-F. Donati 3 and M. Jardine 1 

1 SUPA, School of Physics and Astronomy, University of St Andrews, North Haugh, St Andrews, Fife, KY16 9SS 
2 Department of Astronomy, University of Virginia, P.O. Box 400325, Charlottesville, VA 22904-4325, USA 
3 Laboratoire d'Astrophysique, Observatoire Midi-Pyrenees, 14 Av. E. Belin, F-31400 Toulouse, France 



ABSTRACT 

Models of magnctosphcric accretion on to classical T Tauri stars often assume that 
stellar magnetic fields are simple dipoles. Recently published surface magnetograms of 
BP Tau and V2129 Oph have shown, however, that their fields are more complex. The 
magnetic field of V2129 Oph was found to be predominantly octupolar. For BP Tau the 
magnetic energy was shared mainly between the dipolc and octupole field components, 
with the dipole component being almost four times as strong as that of V2129 Oph. 
From the published surface maps of the photosphcric magnetic fields we extrapolate 
the coronal fields of both stars, and compare the resulting field structures with that of 
a dipole. We consider different models where the disc is truncated at, or well-within, 
the Keplerian corotation radius. We find that although the structure of the surface 
magnetic field is particularly complex for both stars, the geometry of the larger scale 
field, along which accretion is occurring, is somewhat simpler. However, the larger scale 
field is distorted close to the star by the stronger field regions, with the net effect being 
that the fractional open flux through the stellar surface is less than would be expected 
with a dipole magnetic field model. Finally, we estimate the disc truncation radius, 
assuming that this occurs where the magnetic torque from the stellar magnctosphcrc 
is comparable to the viscous torque in the disc. 

Key words: Stars: pre-main sequence - Stars: magnetic fields - Stars: coronae - 
Stars: activity - Xrays: stars - Stars: individual: BP Tau, V2129 Oph 



1 INTRODUCTION 

Classical T Tauri stars (cTTSs) are young solar analogs 
which are accreting material from circumstellar discs. Many 
observations are consistent with the scenario of the stel- 
lar field truncating the inner disc and channelling gas 
onto discrete regions of the stellar surface . The shapes of 
near- IR spectral energy distributions (e.g. iRobitaille et al.l 
120071). and the kinematics of CO lines formed in the disc 
( Naiita. Carr fc Mathieull2003l ) , are consistent with the disc 
having been truncated at a distance of a few stellar radii. 
Average surface fi elds of order a kG h ave been detected on a 
number of cTTSs (| Johns-Krul]||2007f l . which models sugges t 
will be strong enough to disrupt the inner disc (|Koniglfll99lf ) . 
The detection of inverse P-Cygni profiles, with widths of sev- 
eral hundred kms -1 , can also be explained by material es- 
sentially free- falling along the field lines of th e stellar magne- 
tosph ere from the location of the inner disc (JEdwards et al.l 
1 1994 ). Furthermore, the excess continuum emission (veiling) 
in the optical and UV is likely to arise from shocks at the 
base of accretion funnels, with emission lines with high ex- 
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citation potentials, e.g. Hel 5 876A, forming mainly in such 
regions (|Beristain et al.ll2001i ). 

The magnetic interaction between the stellar 
field and the disc may also have important conse- 
quences for the form a tion of planets. Simulations by 
Romanova fc Lovelace! (120061) . and analy tic w o rk b y 



Lin. Bodenheimer fc Richardson! (119961 ) and iFleckl (|200Sh . 



suggest that the inner disc hole, cleared by the star-disc 
interaction, may provide a natural barrier that decreases 
the rate of inward migration of forming planets. There 
is also evidence that the large scale magnetosphere may 
directly disrupt the inner disc, producing warps which 
in some systems c ross the observer's line-of-sight to the 
star (e.g. AA Tau, iBouvier et all 120071 ) . Indeed 3D MHD 
simulations have demonstrated that complicated warping 
effects in th e disc truncation region may be common for T 
Tauri stars (|Romanova et alj|2003t 12004 ). 

The star-disc interaction may also explain the slower 
rotation of cTTSs compared to the typically older weak-line 
T Tauri stars, whose discs h ave largely dispersed (see e.g. 
the review by iBouvieri 120071 ) . Accretion of material from 
the inner disc would act to spin-up the central star in the 
absence of some physical mechanism to remove angular mo- 
mentum from the system. Various magnetospheric accretion 
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models have been developed to explain this mechanism, 
which differ in their assumed location of the inner disc, the 
details of how angular momentum is removed, and how 
the magnet ic field t opolo gy controls both accretion and 
outflows ilKonigl Il99ll; ICollier Cameron fc Campbell 



2000 



1993: IShu et all 1 19941: " iFerreira. Pelletier fc Appj 



Ferrcira, Doug ados fc Cabritl 



2006 



Kiiker Henning fc Rudiger! l2003l: iMatt fc Pudrit3 l2005allbl R c° 



2008al lbl: lLong. Romanova fc Lovelace! 12005 : iBessolaz et al.l 
20081 ). Although observations indicate that the magnetic 



field to pologies of cTTSs are complex (e.g. see the discus- 
sion in iGregorv et al.l l2006aT ). the majority of accretion 
models simply assume that the star has a strong dipolar 
magnetic field. Recently however a few st udies have dr opped 
this assumpt i on, w ith IGregorv et al.l (|2005l . l2006al ) and 
Ijardine et al.l (J2006J) being the first to consider how mag- 
netic fields with a realistic degree of complexity would affect 
the accretion process. Other field geometries hav e since 
been considered by lLong. Romanova fc Lovelace! (J2007I . 
120081 ) who present the results of a 3D MHD simulation 
of accretion to stars with composite fields that consist of 
dipole and quadrupole components tilted at various angles 
relative to the stellar rotation axis. 

It is only recently, however, that instrumentation has 
advanced to the stage where it is possible to map the 
magnetic topo logies of classical T Tauri stars. iDonati et al.l 
42007|, l2008aT ) have recently published magnetic surface 
maps derived from Zeeman-Doppler imaging of the young 
solar-analogs V2129 Oph and BP Tau, using data from 
the ESPaDOnS and NARVAL spectropolarimeters at the 
Canada-France-Hawai'i Telescope and Telescope Bernard 
Lyot respect ively. Using the pre-main seque n ce evolution- 
ary m odel of lSiess. Dufour fc Forestinil (|2000l ). IDonati et al.l 
(|2007h determined that V2129 Oph, a high-mass T Tauri 
star of 1.35 Mo, had developed a radiative core. In con- 
trast, BP Tau, at 0. 7 Mq, should be completely convective 
(JDonati et al.ll2008aO . As the internal structure, and there- 
fore the magnetic field generation process, is different in both 
V2129 Oph and BP Tau, it is of particular interest to com- 
pare the structure and properties of their magnetic fields, as 
well as comparing both to a dipole. 

In §2 we summarise the stellar parameters and mag- 
netic field measurements of BP Tau and V2129 Oph. In §3 
we describe how three-dimensional coronal fields can be ex- 
trapolated from Zeeman-Doppler maps of photospheric mag- 
netic fields, and compare the resulting field topologies with 
a large scale dipole. In §4 and §5 we discuss the structure of 
the open stellar, and accreting, field, while §6 contains our 
conclusions. 



2 THE STARS BP TAU AND V2129 OPH 

BP Tau is one of the best studied classical T Tauri stars, 
with an abundance of observational data across most wave- 
bands. In contrast to this, V2129 Oph is less well studied, 
despite being the brightest T Tauri star in the p-Oph star 
forming cloud. 



2.1 Stellar parameters 

Fo r BP Tau we adopt t he same stellar parameters as used 
bv lDonati et alj (|2008af ) namely a mass of Af* = 0.7 Mq, a 
radius of R„ = 1.95 Rq and a rotation period of P ro t = 7.6 d. 
This implies an equatorial corotation radius, 



GM„P^ 

4tj-2 



1/3 



(1) 



of 7AR*. We assume a mass acc retion rate of 2.88 x 
10- 8 M Q yr- 1 (|Gullbring et al.lll99Sf ). 

Fo r V2129 Oph we ado pt th e same stellar parameters as 
used b v lDonati et al.l (|2007f ) and ljardine. Gregory fc Donatil 
{20081), namely a mass of M* = 1.35 Mq, a radius of 
R* — 2.4 Rq and a rotation period of P ro t = 6.53 d. This im- 
plies an equatorial corotation radius of R co = 6.7 R*. There 
are few estimates of the mas s accretion rate on to V2129 Oph 
available in the literature. Eisner et al.l (|2005T ) . who refer 
to V2129 Oph as AS 207A, estimate an a ccreti on rate of 
3.2 x 10 _8 M Q yr _1 , whereas IDonati et al] (|2007r ) calculate 
a value of 4 x lO^Mgyr" 1 using the empirical relation- 
ship between accr etion rate and the flux in the Cal l 8662A 
line presented by iMohantv. Javawardhana fc Basril l|2005h . 
We therefore assume that the disc mass accretion rate of 
V2129 Oph is A/ = lO^Mgyr -1 , a compromise between 
both observationally inferred values. 



2.2 Magnetic field measurements 



From Zeeman-Doppler imaging of V2129 Oph. IDonati et al.l 
l|2007l ) detected clear circular polarisation signals in both 
photospheric absorption, and accretion related emission, line 
profiles. The temporal variations of the Zeeman signatures 
was dominated by rotational modulation. In deriving a sur- 
face magnetogram it is possible to determine how the mag- 
netic energy is distributed between the different field modes 
of a complex multi-polar field. In other words it is possible to 
determine how strong the dipole component is relative to the 
quadrupole component, or the octupole component, or the 
hexadecapole component etc. In the case of V2129 Oph the 
field is dominated by a 1.2 kG octupole, tilted at ~20°with 
respect to the stellar rotation axis. The dipole component 
was found to be we ak, with a polar str ength of only 0.35 kG 
and tilted at ~30° (|Donati et al.ll2007l '). The surface field in 
the visible hemisphere of V2129 Oph is dominated by a 2 kG 
positive radial field spot at high latitude, with the footpoints 
of the accretion funnel r ooted in this region , but differs sig- 
nificantly from a dipole (JDonati et al.l 120071 ). We note that 
due to the limited resolution achievable by Zeeman-Doppler 
imaging studies, derived surface magnetograms likely miss 
the smallest scale field at the stellar surface. The integrated 
affects of many flares due to magnetic reconnection events 
within such small scale field regions, is likely responsible for 
the "quiesce nt" level of X-ray of emission detected from T 
Tauri stars (|Getman et al.ll2008al ) . 

Several authors have reported the detection of magnetic 
fields on BP Tau dJohns-Krull. Valenti fc Koreskol Il999bl : 
IJohns-KruIil 120071 : IDonati et al.l l2008al ) with a strong cir- 
cular polarisation signal commonly dete cted in the accre- 
tion related Hel 5876A emission line (Jjoh ns-Krul l et al.l 
1999al: IVaTenti fc Johns-Krulll 120041 : ISvmington et al.ll2005l : 
Chuntonov et al.l 120071 ). as well as in other accretion 
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Figure 1. Field extrapolations from observationally derived surface magnetograms of BP Tau (top row) and V2129 Oph (bottom row). 
The left-hand panel shows the complex surface field, the middle panel the more well-ordered larger scale field, and the right-hand panel 
the open field lines. The stellar surface is coloured to show the polarity of the radial field component with red (blue) depicting positive 
(negative) field regions. The images are not to scale. 
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the Call IRT and Hel 
Despite the c oncerns of 
IChuntonov et all (|2007l fl iDonati et all (J2008al ) have pre- 
sented clear evidence that the Zeeman signatures in both 
photospheric and accretion related lines are dominated by 
rotational modulation, having monitored BP Tau for a com- 
plete rotation cycle at two different epochs. Interesting the 
large scale structure of BP Tau's field was found to be sim- 
ilar despite the observing runs being separated by almost 
11 months. However, the magnetic features appear to have 
undergone an apparent phase shift of a quarter of a rota- 
tion phase, which was most likely ca used by a small error in 
the adopted stellar rotation period (|Donati et al.ll2008al )l 2 l 
Although BP Tau's magnetic field is found to be complex, 



1 Some authors have argued that variations in the longitu- 
dinal (line-of-sight) field component, derived from polarisation 
detections in the Hel 5 876A line, were attributable t o ro- 
tational modulation (e.g. | V alenti. Johns-Kr ull fc Hatzed 12003 ; 
IValenti fc Johns-Krunl |2004) . IChuntonov et al.l J2007l ~refutes 

such suggestions and argues that the field in the Hel 5876A line 
formation region is constantly evolving and restructuring on a 
timescale of only a few hours. T he ESPaDOnS/NAR VAL spec- 
tropolarimetric data presented bv lDonati et al.l (I2008ar) . however, 
clearly show that although the Hel 5876A line is subject to in- 
trinsic variability, its temporal evolution is dominated by rota- 
tional modulation. This suggests the magnetic field in the Hel 
line formation region remains stable on timescales of longer than 
a rotation cycle. 

2 Due to the similarity of the de rived surface magneto grams we 
only use the one derived from the lDonati et al,l (|2008aj) February 
2006 data in this paper. 



it is simpler than that of V2129 Oph. The field of BP Tau 
consists of both a strong dipole component, which at 1.2 kG 
is four times as strong as that of V2129 Oph, and a strong 
octupole component of 1.6 kG. Both the octupole and dipole 
moments are tilted by ~ 10° with respect to the stellar rota- 
tion axis, but in different planes (jDonati et al.ll2008al ). 



3 FIELD EXTRAPOLATION 

A variety of field extrapolation techniques have been de- 
veloped which allow the extrapolation of coronal mag- 
netic fields from observationally derived surface magne- 
tograms. From the derived maps of the photospheric fields 
of V2129 Oph and BP Tau, we extrapolate their three- 
dimensional coronal field topologies using the potential field 
source surface model. A potential field is one which is 
current-free, and has the advantages over a full MHD model 
of faster computation speed and simplicity, and does not re- 
quire an assumption about an equation of state. The disad- 
vantages of such a model, however, are that time-dependent 
and non-potential effects cannot be reproduced, and these 
may be important when consi dering the in t eracti on between 
the stellar field and the disc. I Riley et al.l (|2006J ) provide a 
comparison of the two techniques, as applied to the Sun, 
and find that often the potential field mod el produces results 
that closely match MHD models, with iLiu fc Linl (|2007l ) 
coming to the same conclusion for the case of stable so- 
lar active regions. For cTTSs however, the main source of 
non-potentiality in the field is likely to be caused by the in- 
teraction between the stellar magnetosphere and the disc, 
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which cannot be modelled with the static field structures 
considered here. 

Alternative field extrapolation techniques from 
Zeeman-Dopp le r im ages have been considered by 
iHussain et al.l (|2002r ). who compared the difference be- 
tween a non-potential current carrying magnetic field 
model, and the potential field model. By considering 
field extrapolations from a magnetogram of the young 
rapid-rotator AB Dor, they found that the ratio of free 
energy in the potential to that in the non-potential model 
differed by up to 20% close to the stellar surface. However, 
this difference decayed rapidly with height above the star, 
with only 1% difference at 17?* above the photosphere, 
indicating little difference in the field structures derived 
from the potential and non-potential field models. The 
potential field source surface model should therefore be 
adequate to allow us to address the question of whether 
the surface magnetograms of V2129 Oph and BP Tau are 
consistent with the observed locations of accretion funnels 
and accretion hotspots, as well as to compare the large 
scale field topologies with that of a simple dipole. 

As we have assumed that the field is potential then 
V x B — 0. This condition can be satisfied by writing the 
field in terms of a scalar flux function *]/ 



-V*. 



(2) 



The field must also satisfy Maxwell's equation, V • B = 0, 
and therefore "$> must satisfy Laplace's equation, V 2 * = 0. 
This has a standard solution of the form, 



* = J2 J2 [ aimr ' + hmr -C,+1) ] Pim(e)e™ 



(3) 



Z = l m — — I 



where Pi m denote the associated Legendre functions and the 
coefficients ai m and bi m are determined from the boundary 
conditions. The first boundary condition is that the strength 
of the radial component of the field at the stellar surface 
is that which is derived from the Zeeman-Doppler maps 
of the star. The second condition is that at some height 
above the stellar surface Rs, known as the source surface, 
the field becomes radial and hence Be(Rs) = B$(Rs) — 0, 
emulating the effect of the corona blowing open fie ld lines 
to form a stellar wind (jAltschuler fc NewkirkHl969T ). Thus, 
from ([2]) it is possible to determine the magnetic field compo- 
nents, B r , B$, and B$, and therefo re a field vector , at any 
point within the stars corona (see iGregorv et al.l (|2006al ) 
for the full expressions). In order to extrapolate the field 
we used a modified version of a code origina lly developed 
bv lvan Ballegooiien. Cartledge fc Priest (1199a). and ap plied 
specifically to T Tau ri stars by Ijardine et al.l (|2006l ) and 
IGregorv et all i|2006al lbl. 120071 ) . 

The field extrapolations showing the smaller scale, 
larger scale, and open field structure of V2129 Oph and 
BP Tau, are shown in Fig. [T] For both stars we have set 
the source surface to be at the corotation radius (see §2.1), 
which at this stage serves only to illustrate the difference be- 
tween the complex and loopy surface field regions, and the 
well-ordered larger scale field. Later we discuss how varia- 
tions in the source surface radius affects the structure of the 
field. It is clear from Fig. Q] that for both stars the largest 
scale field resembles a slightly tilted dipole, especially in the 
case BP Tau. But how does the larger scale field extrap- 



olated from observationally derived surface magnetograms, 
compare to a simple dipole magnetic field? 

3.1 Comparison with a large scale dipole field 

The complex field structures obtained by field extrapolation 
of both V2129 Oph, and in particular BP Tau, appear to 
indicate that the larger scale field is well ordered, and much 
simpler than the loopy and compact surface field regions. 
In this section we compare the larger scale field topology 
with that of a modified dipolar field. Taking the I — 1, m — 
component of equation © and imposing the boundary 
conditions 



B r (R„) = |^cos(9 



B e (Rs 



B^Rs 



0. 



(4) 
(5) 



allows the construction of a dipole field with a source surface, 
where fi — i? H ,_B* :PO i e /2 is the dipole moment and B t:PO i e 
the polar strength of the dipole. The inclusion of a source 
surface introduces a modification to the standard dipole field 
components, 



B r — 

Be = 



2^ cos 9 ( r 3 + 2R% 



Rl + 2R% 



/ism( 



0. 



-2r 3 + 2R 



which are recovered in the limit of Rs 
lines can be described by, 

sin 6 — — — =-, 

r 3 + 2R% ' 



(6) 

(7) 

(8) 
oo. Individual field 



(9) 



where £ is a constant along a particular field line, such that 
different values of £ correspond to different field lines. The 
last closed field line passes through the stars equatorial plane 
(6 = 7r/2) at r = r m = Rs, and therefore £ = 3R%. Such a 
field line connects to the stellar surface at a co-latitude of 
8 m where 



dRcjR* 

rJT2rJ 



(10) 



The struct ure of such a modified dipolar field is illustrated 
in fig. 1 of Ijardine et al.l (|2006r ). We note that the dipole 
components of the fields of both BP Tau and V2129 Oph are 
slightly tilted with respect to the stellar rotation axis (10° for 
BP Tau and 30°for V2129 Oph). We account for this small 
tilt below when comparin g the complex field s t ructu res to 
that of the dipolar fields. iMahdavi fc Kenvonl (|1998l ) pro- 
vide analytic expressions for calculating the components of 
an arbitrarily tilted dipole, which we adapt in order to in- 
clude the source surface, and calculate the tilted dipole field 
components numerically using our field extrapolation code. 
The black points in Figs. [2] and [3] show field line length 
against width calculated from the field extrapolations of the 
magnetic fields of BP Tau and V2129 Oph, for source surface 
radii of Rs = 2.5 R, and Rs ~ R co - At this stage such plots 
are used only to illustrate by how much the field structure 
departs from a simple dipole (red points in Figs.[2]and|3J). We 
define the length of a field line as being the radial distance 
measured from the stellar surface to the maximum height 



The non-dipolar fields of T Tauri stars 5 



1.5 



a l 



c 
v 

0.5 





II 1 1 1 1 1 II 


i i 


III 




BP Tau 






- 


R S «2.5R„ |j 

1 

■ ■ 

1 

■ i 


■ : ^'v". i. 


1 ' : >7*& 




! 1 

■ ■ 

■ I 




*w 


- 


/ 




lip? 


- 


■ i 
•1 A 


M&i.-'.'-f-. : 


jjtiS 




i i i i i i i i i 


I I I 


1 1 1 1 1 1 



1.5 



03 1 



ao 

c 

0.5 



1 


1 


1 1 1 1 1 


iii 


i i i i 







V2129 Oph , 









- 




J 


ipftK: 




■ 


- 


R s *2 5K 


t ■ 




- 














- 






1®^^^'" * 




- 


- 






' 




- 


- 




4 

1 

la 


P^ . 


if;';-'' ^' c -' 


- 


- 






P^^Sp^"' ; J 


4^' ■'. 


- 






$M 


|ajF r / 










Bffl 












■■jsp 








- 




0Jgm$' 






- 




l l 


I i i i i 


1 i i i 


, i , , 


i i 



0.5 1 1.5 

width/R + 



0.5 1 1.5 

width/R + 



Figure 2. Plots of field line length (maximum radial distance above the stellar surface) against field line width (the distance along a 
segment of the great circle connecting the field line footpoints on the stellar surface) for BP Tau (left panel) and V2129 Oph (right 
panel). A source surface radius of 2.5 iJ* has been assumed. The black points are for the extrapolated fields, the red points indicate the 
behaviour of a dipole field, and the green points an octupolc field. 
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Figure 3. The same as Fig. [2] but assuming that the source surface is at approximately the corotation radius, which is 7.4 R* for BP Tau 
(left panel) and 6.7 R* for V2129 Oph (right panel). 



of the field line. The width of a field line is defined as the 
distance along the segment of the great circle connecting the 
field line footpoints on the stellar surface. For a field line 
with footpoints at co-latitudes and longitudes (8± , (pi ) and 
(82, 02) the width, w, can be calculated using the Haversine 
formula. 



h = 



I I + sm oi sm 02 sm I I 



V 2 J 



V 2 ) 



(11) 



= 2sin _1 (\/ft) 



(12) 



where A<j) — (f>2 — <p\ is the longitudinal separation between 
the footpoints, A9 — |#2 — #i| is the difference in co-latitudes, 
and where all angles are measured in radians and all dis- 
tances in units of stellar radii. 

By considering progressively longer field line loops, Figs. 
[2] and [3] show that for a dipole field the field line width 
eventually tends asymptotically to a certain value which de- 
pends on the choice of source surface. This is as expected, 
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as by equation (|10[) the co-latitude of the footpoint of the 
longest dipole field line depends on the choice of source sur- 
face Rs, with 6 m being smaller for a larger Rs (and for an 
aligned dipole field line the latitudinal footpoint separation 
is A9 — tt — 29m radians). We also show in Figs. [2] and [3] 
the behaviour of an octupole magnetic field (green points). 
For an axial octupole (I — 3, m = 0) the width of the closed 
field line loops about the stellar equa tor is less than those a t 
higher latitudes (see the discussion in lWillis fc Younal 19871 ). 
which gives rise to two separate lines in the length vs width 
plots. 

For a small source surface (see Fig. the longer 
field lines on BP Tau are typically wider than those on 
V2129 Oph. This reflects the strength of the dipole com- 
ponents in the two stars, which is about 4-times stronger on 
BP Tau than V2129 Oph (stars with a strong dipolar field 
component would typically have wider field lines, connect- 
ing the different polarity poles in opposite hemispheres), and 
also the dominantly octupolar nature of V2129 Oph's field. 
In both cases, however, the fields are more complex than a 
dipole. For the case of a large source surface (see Fig. |3J) the 
field structure of BP Tau follows a similar trend to the dipole 
field, as do the longest field lines on V2129 Oph. However, 
many of the smaller scale field lines on V2129 Oph are not as 
wide as the dipolar field lines. This suggests that the smaller 
scale field on V2129 Oph is complex and multi-polar with 
the field strength decaying more rapidly with height than on 
BP Tau; we return to this point in §5. For V2129 Oph the 
smaller scale field (see Fig. [2j is more closely matched to an 
octupole field. Also, some of the shortest width field lines 
are longer than those of an octupole, indicating a contribu- 
tion to the surface field from the higher order field modes. 
Indeed, for V2129 Oph there is also a s ignificant amount o f 
magnetic energy in the I = 5 field mode (|Donati et al.ll2007h . 

The field structure of V2129 Oph is particularly strik- 
ing when a large source surface is considered (Fig. [3] right 
panel). The surface magnetic field of V2129 Oph is domi- 
nated by large 2 kG positive radial field spot at high latitude 
in the visible hemisphere, and li kely a similar negat ive field 
region in the unseen hemisphere (|Donati et al.ll2007l ) . While 
the surface field region consists of many small width loops, 
the larger scale field comprises of the longer magnetic loops 
which connect to the high latitude field regions. This be- 
haviour is not apparent when considering V2129 Oph with 
a smaller source surface, as the co-latitude at which the last 
closed field line connects to the stellar surface is larger, by 
equation (|10|) . which means that field lines cannot connect 
to the large high latitude radial field spots. With a large 
source surface the longest field lines are always wider than 
those of a dipole, suggesting that the structure of the field 
lines are being affected by the strong field regions close to 
the stellar surface; with such field lines connecting at smaller 
co-latitudes than a dipole field. This effect is more appar- 
ent for V2129 Oph, where the larger scale field lines are 
"squeezed" by the strong dominantly octupolar field close 
to the stellar surface. This forces such field lines to have 
larger widths at a given length. This is also the case for the 
field of BP Tau (see Fig|3j left panel), however, the effect is 
less due to the strength of octupole component relative to 
the dipole component being over 2.5 times less for BP Tau 
than for V2129 Oph. The structure of the largest scale field 
lines therefore depart from the path followed by purely dipo- 
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Figure 4. The variation of the ratio of open to total flux through 
the stellar surface for a range of source surface radii. The solid 
line represents the modified dipole magnetic field, the dashed line 
the field of BP Tau and the dotted line the field of V2129 Oph. 
There is less fractional open flux when considering magnetic fields 
with a realistic degree of complexity. 



lar field lines. It is such larger scale loops which are likely 
to be carrying material from the disc to the star. Therefore 
the location of accretion hotspots on the stellar surface will 
differ from what is expected from dipolar accretion, even 
though the larger scale field is relatively well-ordered and 
similar to a dipole. In the following sections we look more 
closely at the structure of the open and accreting field. 



4 THE OPEN FIELD 

The amount of open flux passing through the stellar surface 
is important for models of stellar winds, which for cTTSs 
may remove the angular angular moment um which would 
be tra nsfered to the star due to accretion (|Matt fc Pudrita 
l2005br i. It is therefore of interest to compare the amount 
of open flux for V2129 Oph and BP Tau with that of the 
modified dipole magnetic field. Since the magnetic torque in 
a stellar wind does not depend on the polarity of the field, 
we consider here the absolute value of the flux. The open 
flux through the stellar surface is given by 



\B r , open (0,(/))\dQ. 



(13) 



Using equation (4) an expression for $ opE „ can be derived 
for our modified dipole field, 



2[i 
RS 



— 2irR t B t ,„ ; e sin 



(14) 
(15) 



where the additional factor of 2 accounts for the open field 
in each hemisphere. A similar expression can be derived for 
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Figure 5. The variation of open and closed flux (left and right panels respectively) through the stellar surface for a range of source 
surface radii. The solid lines represent modified dipole magnetic fields, with the red (black) line representing a polar field strength of 1.2 
(0.35) kG. The dashed line represents the field of BP Tau and the dotted line the field of V2129 Oph. There is less open flux, but more 
closed flux, when considering the complex magnetic fields. 



the closed flux through the stellar surface, 



^closed 



■1R 



/ '* „, T /2 2/i 



do 



m 



cos 9 sin 6d6 



2-kR 2 ,B,, p l e COS Or, 



(16) 

(17) 



Thus for the modified dipole field the ratio of open flux to 
the total flux through the stellar surface is, 



<l>o 



$< 



3x^5 XI* 

m + 2R 3 s ' 



(18) 



where the total flux $ to tai = ®o P en + ^closed- The solid 
line in Fig. [4] shows the variation of the fractional open flux 
&o P en/&totai for a range of different source surface radii Rs 
assuming a modified dipolar magnetic field. For large source 
surface radii the colatitude of the footpoint of the last closed 
field line is close to the pole of the star, and therefore the 
open flux through the stellar surface is less. As the source 
surface radius is decreased more of the closed field is con- 
verted to open field along which a (stellar) wind may be 
launched. 

For the magnetic fields extrapolated from the surface 
magnetograms we divide the stellar surface into a series of 
small grid cells. Within each grid cell we take the average of 
the magnitude of B r , at the footpoint of the open field lines 
and multiply this by the area of the cell in order to calculate 
the open flux from that individual cell. If a cell contains the 
footpoints of both open and closed field lines we weight the 
area of that cell accordingly. For example, a cell containing 
50% open field, is assumed to contribute 50% of it's area 
to the open flux. The total open flux through the stellar 
surface is then obtained by summing over all the cells, while 
the total closed flux is calculated in a similar way. We find 
that there is less fractional open flux through the surface of 
V2129 Oph and BP Tau (dotted and dashed line in Fig. 0} 



compared with the dipole field. This is consistent with our 
argument in the previous section that the structure of the 
larger scale field (and in particular that of V2129 Oph) is 
influenced by the strong field regions closer to the star. As 
the larger scale field lines are wider than those of a dipole, 
their footpoints are closer to the pole (see Fig. [3}. Thus 
there is less area of the stellar surface available from which 
a stellar wind can be launched along the open field. 



The lower fractional open flux for the complex fields, 
can be understood by considering how the open and closed 
flux changes as the source surface location is varied. This is 
illustrated in Fig. [S] where the solid lines represent dipole 
magnetic fields. In order to compare the fluxes of the com- 
plex fields with a dipole, we have assumed that the dipole 
has a polar strength equal to that of the measured dipole 
components of V2129 Oph and BP Tau (0.35 kG and 1.2 kG 
respectively). Both stars have slightly less open flux com- 
pared to a dipole, but more closed flux. There is significantly 
more closed flux for V2129 Oph compared to a dipole than 
for BP Tau. This is simply a reflection of the more complex 
nature of V2129 Oph's magnetic field, and is the largest fac- 
tor in explaining why the fractional open flux for this star 
is well below the dipole case (see Fig. \%§. 



The main result of this section is that there is less frac- 
tional open flux for the complex field geometries, compared 
to a pure dipole case. This is result of the complex fields hav- 
ing both lower open flux, and, in particular, higher closed 
flux. This is understandable since a more complex field 
geometry is characterised by more numerous and smaller 
closed magnetic loops. 
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5 THE ACCRETING FIELD 

It has long been suspected that the magnetic fields of cTTSs 
are more comp l ex th an a simple dipole, with, for example, 
lEdwards et al.l (|1994f ) arguing that the common detection 
of inverse P-Cygni profiles was evidence for multiple ac- 
creting loops in a field configuration more complex than a 
dipole. Until very recently, however, only indirect clues were 
available to the complex nature of T Tauri magnetic fields. 
The detection of rotationally modulated X-ray emission in 
a small, but significant, sample of T Tauri stars indicated 
that X-ray emitting plasma was confined to compact field 
regio ns in-homogeneously d istributed about the stellar sur- 
face (|Flaccomio et alj2005f ) . Such modulated X-ray emission 
would not occur if the surface field was dipolar. Further clues 
to the complex nature of T Tauri surface fields comes from 
the common failure to detect polarisation signals in pho- 
tospheric absorption lines, despite strong average fields be- 
ing measured from the Zeeman broadened, and unpolarised, 
intensity spectra (jValenti fc Johns-Krulll 12004 ). This sug- 
gested that the surface of T Tauri stars were covered in many 
regions of opposite polarity. Spectral lines forming in such 
magnetic regions would be polarised in the opposite sense, 
leading to a cancellation of the polarisation signal along the 
line-of-sight to the observer. However, often a strong polar- 
isati on signal is detected in a ccretion related emission lines 
fe.g. ljohns-Krull et alll999ah which is modu lated due to the 
stellar rotation (|Valenti fc Johns-Krullll2004r ). This suggests 
that such emission lines are forming in regions of the stellar 
atmosphere which are permeated by field lines of a single 
magnetic polarity which remains stable as the star rotates. 
In other words, despite the clues to the complex nature of 
T Tauri surface fields, the accreting field lines are likely an- 
chored in single polarity regions, and the field connecting to 
such regions may be well-ordered and more like a dipole. 

The resolution achievable with ESPaDOnS (and its twin 
instrument NARVAL) has allowed the detection of polari- 
sation signals in both photospheric absorption, and accre- 
tion related emission, lines, and has c onfirmed the com- 
plex n ature of magnetic fields of cTTSs (|Donati et alj|2007l 
l2008aT ). However, magnetospheric accretion models devel- 
oped thus-far have yet to account for the true complexity 
of the stellar magnetic field. The first steps in consider- 
ing how realistic multi-polar fields w ould influence accre- 
tion fr om t he inner disc were ta ken bv lGregorv et al.l (|2005l . 
1200681 ) and lJardine et alj (120061). More r ecently the 3D MHD 
simulations of lLong et al.l ( 20071 . 120080 have considered a 
quadrupole-dipole composite field. However, simulations in- 
corporating generalised multi-polar magnetic fields are now 
required, and will shed new-light on the star-disc interaction. 

For the stars where magnetic surface maps have been 
obtained to date the bulk of accretion in the visible hemi- 
sphere i£_cliarnTfiU^ im^a smgle rnagnetic spot at high lat- 
itude (|Donati et all 120071 . l2008ah . In order to explain the 
location of high latitude hotspots, the source surface has to 
be large enough to allow field lines to connect between the 
high latitude radial field spots close to the poles. The source 
surface radius represents the extent of the closed stellar mag- 
netosphere and independently obtained observations provide 
justification for setting t he source surface to be at least equal 
to the corotation radius. iGetman et al.1 ( 2008a. b) have anal- 
ysed flares detected during the Chandra Orion Ultradeep 



Project, using a new technique that allows the modelling of 
flaring events which are much fainter than those accessible 
by traditional par ametric flare modelling techniques (e.g. 
iFavata et al.ll2005r ). By calculat ing the lengths of the mag- 
netic loops containing the flares. IGetman et al.1 (|2008br ) have 
found evidence for magnetic structures that extend well be- 
yond the corotation radius radius in weak-line T Tauri stars, 
however for classical T Tauri stars the magnetic loops are all 
confined to within the corotation radius. This suggests the 
presence of a disc is limiting the extent of the closed stellar 
magnetos phere of accret i ng T T auri s tars - as previously pre - 
dicted by iJardine et al.1 (|2006l ). The IGetman et al.1 (|2008bl ) 
results therefore provide strong evidence that the stellar 
magnetosphere in classical T Tauri stars never extends to 
beyond the corotation radius, suggesting that setting the 
source surface to be at corotation is a reasonable assump- 
tion. However, the location of the disc truncation radius may 
still be closer to the stellar surface. With a large source sur- 
face field lines are able to connect to the high latitude field 
regions, even well within the corotation radius. Thus if the 
disc is truncated at, or even within, the corotation radius, 
then accretion may proceed along field lines which connect 
to a high enough latitude in order to explain the observed 
hotspot locations. Therefore an interesting question to ask is 
how would the disc truncation radius differ between dipole 
and more complex magnetic field models? 

5.1 The inner disc radius 

The question of where the disc is truncated, and therefore 
the structure of the magnetic field threading the disc 
and carrying accreting gas, remains a major problem 
for accretion models. It is still unknown if the disc is 
truncated in the vicinity of the corotation radi us, the 
assum ption of traditional accretion models (e.g. iKonigll 
Il99ll), or whether it extend s closer to the stellar s urface 
(e.g. iMatt fcPud"rit3l2005bl ). As in iGregorv et all (|2007f ) 
we assume that accretion occurs over a range of radii 
within the corotation radius . This is equivalent to the 
appro ach tak en previously by iHartmann. Hewett fc Calveti 



Il994l). iMuzerolle. Calyet fc Hartmannl J20011 ). 

Symington et al.1 (120051). lAzevedo et al.1 (|2006T ) and 



Kurosawa. Harries fc Symington! ( 20061 ) who have demon- 
strated that such an assumption broadly reproduces 
observed spectral line profiles and variability. It should 
also be noted that the accreting field geometries which we 
consider here are only snap-shots in time, and in reality will 
evol ve due to the interact ion with the disc. 

iBessolaz et al.1 (|2008l ) provide an overview of the var- 
ious assumptions that have been used in the literature to 
calculate the disc truncation radius. Here we provide rough 
estimates of the inner disc radius for BP Tau and V2129 Oph 
based on their extrapolated magnetic fields, and compare 
these values with those calculated by assuming both stars 
have dipole magnetic fields. 

We assume that the location of the inner disc is the ra- 
dius at which the torque due to viscous processes in the disc 
is comparable to the magnetic torque due to the stellar mag- 
netosphere. This is equivalent to the a pproach taken previ - 
ous ly by several authors, for example, IClarke et al.1 (J1995T ) 
and lWand (|l996T ). This torque depends on how much the disc 
is able to twist the field lines of the stellar magnetosphere. 
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Figure 6. The variation along the equatorial plane of the differential magnetic/viscous torque [i.e. the LHS and RHS of equation | |19| |1 
for BP Tau (left panel) and V2129 Oph (right panel). The solid line represents the RHS of 1190 . with the LHS for an aligned dipolc field 
(dashed line), a tilted dipole field with a source surface (dash-dot line) and the field extrapolation of V2129 Oph/BP Tau (dotted line). 
In both cases the disc is truncated well within the corotation radius of 6.7 R* (V2129 Oph) and 7.4 ij* (BP Tau), although the exact 
location of Rt is very sensitive to the assumed accretion rate. 



We approximate that this twist at the location of the inner 
disc is of order 45° i.e. at the inner disc we assume that the 
perturbed toroidal component is equal to the poloidal com- 
ponent (B<f, ~ B z ). This assumption should be valid as long 
as the disc is truncated sufficiently within the corotation ra- 
dius and assu ming the disc is strongly coupled to the stellar 
magnetic field (|Matt fc Pudritdl2005ah . Equating the differ- 
ential magnetic and viscous torques, and assuming that the 
poloidal field threading the disc is dominated by the vertical 
component (B z 3> B r ), gives 



i*( 



GhU \ 1/2 



(19) 



Note that this equation assumes a cylindrical corordinate 
system, however as we are restricting the problem to the 
star's equatorial plane, the cylindrical radius and the spher- 
ical radius are the same, and B z = Bg. Interestingly, the 
disc disruption ra dius found in the numerical simulations of 
lLong et al.l l|2005l ) (see the discussion in their section 2) ap- 
proximately conincides with that derived through the use of 
equation (|I9|) . The solid line in Fig. [6] represents the varia- 
tion in the RHS of this equation along the stars equatorial 
plane, using the stellar parameters given in §2.1. The various 
dashed and dotted lines in Fig. [6] represent the LHS of this 
equation assuming different forms for the stellar magnetic 
field. The value of r where the lines representing the LHS 
and RHS of equation (|19|) cross, is then the disc truncation 
radius Rt. 

For a dipole magnetic field the variation in B z along 
the equatorial plane can be written as 



B z =B e = -B*, poU [-y i 



(20) 



The dashed lines in Fig. [6] represent such an aligned dipole 



field with a polar strength of B*, po (e = 1.2 kG for BP Tau 
and 0.35 kG for V2129 Oph, i.e. dipole fields with the same 
polar strength as the dipole components of the measured 
fields of both stars (see §2.2). With such an assumption for 
the stellar magnetic fields, we obtain a disc truncation ra- 
dius of ~ 5R* for BP Tau and ~ 3.6R, for V2129 Oph. 
This is well within the corotation radius in both cases, jus- 
tifying our assumption that the perturbed toroidal compo- 
nent of the field equals the poliodal component at the inner 
edge of the disc (Rt ~ 0.7R CO for BP Tau and R t ~ 0.5R CO 
for V2129 Oph). However, for a better comparison with the 
extrapolated fields we consider modified and tilted dipole 
fields, i.e. dipole fields with a source surface, tilted by the 
same amount as the dipole component in both stars. Once 
the dipole field been tilted, the magnetic field is no longer 
axisymmetric in the stars equatorial plane, and therefore 
B z is no longer uniform in azimuth at a fixed radius. At 
each radius r we therefore take the average of the square of 
B z for use in equation (|19[I 3 I (note that we follow the same 
procedure for the extrapolated fields discussed below). The 
dash-dot lines in Fig. [6] represent such tilted dipoles, with 
the disc truncation radius being slightly closer to the stellar 
surface for such fields due to the non-uniformity of B z in 
azimuth. In Fig. [6] (and Figs. [8] and [7] below) we have chosen 
to set Rs = 20 R, in order to minimise the effects of the 
source surface boundary condition upon the field structure 
close to the accreting field regions. 

We have already seen in §3.1 that the larger scale field of 
BP Tau is similar to a dipole, but that the field of V2129 Oph 



3 Once the dipole has been tilted the field threading the disc 
has a small radial component, however for small dipole tilts this 
remains much smaller than B z . 
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shows significant departures from such a simpie field config- 
uration. We therefore expect that the disc truncation radius 
calculated using the extrapolated field of BP Tau will be 
similar to that of a dipole while Rt for V2129 Oph may be 
different. The dotted lines in Fig. [6] demonstrate that this is 
indeed the case, with Rt for BP Tau closely matched to that 
of the modified dipole. This is due to the strong dipole com- 
ponent of BP Tau's field, which c ontains 50% of the mag- 
netic energy (|Donati et alJl2008aD . and is 4-times stronger 
than the dipole component of V2129 Oph. The calculated in- 
ner disc radius for V2129 Oph is closer to the stellar surface 
than would be expected for a dipole field. This is a reflection 
of the intrinsic field complexity of V2129 Oph, and the rapid 
drop off in field strength with height above the stellar surface 
for its magnetic field (see §5.2 below). When the magnetic 
field of the central star is particularly complex, as in the 
case of V2129 Oph, we find that the inner disc is truncated 
closer to the stellar surface than would be expected from 
dipole magnetic field models. However, the use of equations 
such as (|19|) in calculating disc truncation radii for individual 
stars should be treated with caution. Firstly, the values of 
Rt determined here are very sensitive to the assumed mass 
accretion rates M, which are poorly constrained observa- 
tionally. For example, only a factor of 2 difference in our 
assumed accretion rate for BP Tau of 2.88 x 10 -8 M0yr _1 
would change the inner disc radius from 4.87?, to 3.97?* or 
5.97?, depending on whether the accretion rate is increased 
or decreased respectively. Secondly, equation (fl9|) assumes 
that the field threading the disc is dominated by the vertical 
component B z , and although this is the case for the dipole- 
like large scale field of BP Tau, it is not true for V2129 Oph 
where the field has a significant radial component at the 
stellar equatorial plane. Thus the shear within the disc of 
V2129 Oph may generate a significant toroidal component, 
invalidating the use of equation (|19|l . and suggesting that 
our qualitative calculation of the disc truncation radius is 
too simplistic. An MHD simulation using the extrapolated 
fields as starting points would be welcome here, in order 
to calculate inner disc radii more quantitatively. However, 
our conclusion that the inner disc should be truncated at, 
or within, the location obtained when considered a dipole 
magnetic field should remain valid. 

5.2 Structure of the accreting field 

In the previous sections we compared the structure of the 
magnetic fields of V2129 Oph and BP Tau to that of a dipole. 
In this section we select out only those field lines which 
would be able to interact with the accretion disc and carry 
material on to the star. In order to select such fie l d lines 
we fo llow the algorithm discussed bv lGregorv et al.l (|2006al , 
120071 ) . which we only summarise here. 

We consider a thin accretion disc and for V2129 Oph 
assume that the disc normal is parallel to the stellar rota- 
tion axis, such that the disc mid-plane is aligned with the 
stars equator i al plane . By extrapolating the field of BP Tau, 
iDonati et al.l (J2008a|) argued that a flat disc model could 
not explain the observed location of accretion hotspots. 
For material to accrete into the observed hotspot loca- 
tions, either the disc must be tilted relative to the stellar 
rotation axis (i.e. aligned with the stellar magnetic equa- 
tor - similar to the scenario found during the 3D MHD 



simulations of iRomanova et al.l 120031 ) . or the inner disc 
must be warped. Radio CO maps of BP Tau's disc sug- 
gest a disc i nclination of 30°offset from the st ellar incli- 
nation of 45° (|Simon. Dutrev fc Guilloteaull2000l '). although 
both values are rather uncertain. Analytic arguments sug- 
gest that a large scale tilted magnetosphere (as possessed 
by BP Tau) will warp the inner disc, allowing material 
to accrete on to latitudes which are inaccessible if ac- 
cretion were to proceed from a flat, wedge-shaped, disc 
JTerquem fc Papaloizoul l200Ct). Furth ermore, the 3D MHD 
simulations of lRomanova et al.l (|2003T ) , suggest that compli- 
cated inner disc warps may arise due to the interaction of the 
magnetosphere with the disc. There is also strong evidence 
for the stellar field warping the disc in at least one other 
star, AA Tau (e.g. iBouvier et al] 120071 ). which is inclined 
such that our line-of-sight to the star looks through the in- 
ner disc. Given the smaller inclination of BP Tau (i = 45°), 
however, photometric eclipses due to a dense inner disc warp 
crossing our line-of-sight are not expected. For BP Tau we 
therefore assume that the disc is slightly tilted, or equiva- 
lently that the inner disc has been warped by the stellar mag- 
netospher jj. An obv ious question to ask, however, is why did 
IDonati et al] (|2007f ) find that a flat disc model was sufficient 
to explain the observ ed hotspot locations on V2129 Oph, but 
IDonati et al] (|2008aT ) find that a small inner disc warp was 
required for BP Tau? This is likely related to the strength of 
the dipole component on both stars, which is almost 4-times 
stronger on BP Tau than V2129 Oph. At the inner disc, at a 
distance of a few stellar radii, the dipole component is likely 
dominant, with the strength of the higher order field compo- 
nents dropping faster with height above the star. The larger 
scale magnetosphere of BP Tau is therefore more likely to 
distort the circumstellar disc. Indeed to explain the hotspot 
locations on BP Tau the disc had to be tilted in the same 
direction as the dipole moment. Tilting the disc in a similar 
way for V2129 Oph makes little difference to the hotspot 
location, suggesting that it is the strength of the dipole that 
is affecting the inner disc. Further numerical simulations of 
how magnetic fields with a realistically complexity can dis- 
tort the structure of the disc will be required in future to 
confirm this, and are currently being undertaken. 

In order to compare the structure of the complex fields 
of V2129 Oph and BP Tau with that of a dipole we consider 
how fast the field strength drops with height above the stel- 
lar surface. We therefore compare how the ratio B(r)/B, 
varies with r, where B (r) is the strength of the field at some 
position r along the path of a field line and 73* the field 
strength at the field line footpoint. r is the spherical radius 
measured from the centre of the star. Figs.[7]and|S]show such 
plots for both V2129 Oph and BP Tau. All B(r)/B* values 
are plotted from the disc mid-plane to the field line footpoint 
at the stellar surface. Accretion is assumed to occur along all 
field lines threading the mid-plane of the disc across a range 
of radii. In Fig. [7] we assume that accretion is occuring from 
the corotation radius, down to 0.75 7? co (~ 5 — 6.7 7?, for 
V2129 Oph and ~ 5.5 - 7.4 7?, for BP Tau), whereas in 
Fig. [H] we assume that accretion is occuring from a range of 



4 We note, however, that the differences in the accreting field 
structure shown in Figs.[8]and[7]for BP Tau are almost negligible 
when considering a flat or slightly warped disc. 
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Figure 7. The drop in field strength with height above the surface of BP Tau (left panel) and V2129 Oph (right panel) along the 
accreting field lines (black points). Also shown for comparison is a dipole field (red points). The dipole has been tilted by 10° (30°) for 
comparison with the field of BP Tau (V2129 Oph). Accretion is assumed to occur from the corotation radius to ~ 0.75 R co - 
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Figure 8. The same as Fig. [jj but assuming that accretion is occuring from a range of radii about the disc truncation radius, as 
calculated in §5.1 (see Fig. [6}. The structure of the accreting field is shown for BP Tau (left hand panel) and V2129 Oph (right hand 
panel). 



radii about the disc truncation radius calculated in §5.1. For 
V2129 Oph, where the disc truncation radius was calculated 
to be 2.8 Rt, (see Fig. |SJ), the accreting field is assumed to 
span the range of radii ~ 2.4 — 3.2 i?*; for BP Tau, with a 
disc truncation radius of ~ 4.8 R,, accretion is assumed to 
occur from ~ 4.4—5.2 R t . In practise little difference is found 
in the structure of the accreting field when these values are 
varied. Figs. [7J and [5] also show the behaviour of a dipole 
magnetic field, tilted by 10° for comparison with the BP Tau 



data (whose dipole component is tilted by 10°relative to the 
stellar rotation axis), and by 30° for comparison with the 
V2129 Oph data. 

Plots such as Figs. [7] and [S] allow a comparison of how 
the field complexity varies with height above the stellar sur- 
face. We note, however, that we do not account for how the 
structure of the magnetic field will evolve in time due to the 
interaction with the disc. However, this will be less of an is- 
sue if accretion proceeds from close to corotation, where the 
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differential rotation of field line footpoints anchored in the 
disc and on the stellar surface is much less than if the disc 
is truncated well within corotation. As the strength of the 
dipole component decreases with height more slowly than 
the higher order field components, the behaviour of the field 
on the largest scales should be somewhat simpler than that 
of the surface field. This is the case for the magnetic field of 
V2129 Oph where the field behaves more like a dipole above 
~ 2.57?* (see Fig. [7] right panel), after an initial fast drop 
within the strong complex field region close to the stellar 
surface. For BP Tau the field is generally more dipolar, even 
closer to the stellar surface (see Fig. [7J left panel) , reflect- 
ing the much stronger dipole component of it's magnetic 
field, found to b e 1.2 kG and contain ing about 50% of the 
magnetic energy (|Donati et al.ll2008ah . For V2129 Oph if we 
assume that accretion occurs from around the disc trunca- 
tion radius as calculated in the previous section, see Fig. [8] 
we find that the structure of the ac creting field i s much 
more complex than a simple dipole. iDonati et al.l (|2007l ) 
found that the bulk of accretion in the visible hemisphere of 
V2129 Oph occured into a single mag netic spot at high lati - 
tud e. Using the algorithm discussed in lGregorv et al.l (|2007r ) 
and lJardine et al.l (|2008f ) we calculate that about 45% of the 
mass supplied by the disc would accrete into the high lat- 
itude northern hemisphere spot on V2129 Oph, with a si- 
miliar amount to an unseen southern hemisphere spot, and 
less than 10% to smaller equatorial spots. Thus, provided 
that the source surface is set to a large enough radius to 
allow the longest field lines to connect to the high latitude 
spot, the bulk of the accreting material accretes into the 
high latitude spots, even if the disc is truncated within the 
corotation radius. 



6 CONCLUSIONS AND DISCUSSION 

By extrapolating the fields of two T Tauri stars from ob- 
servationally derived surface magnetograms we have com- 
pared the resulting magnetic field topologies with a sim- 
ple dipole. We have found that although the surface mag- 
netic fields of both stars are particularly complex, the larger 
scale field is simpler and more well ordered. This is consis- 
tent with previously published spectropolarimetric studies 
of accreting T Tauri stars, which suggested that bundles 
of accreting field lines connect to the stellar surface in sin- 
gle polarity regions, and were th erefore likely well-ordered 
(e.g. IValenti fc Johns-Krullll2004 l. The larger scale field of 
BP Tau in particular is closely matched to a dipole mag- 
netic field, although V2129 Oph still shows departures from 
a dipole field even on the largest scales. However, for both 
stars, the footpoint separation at the stellar surface of the 
largest scale field lines is greater than for a purely dipolar 
magnetic field. In other words as the largest scale field lines 
arrive at the stellar surface, their structure is influenced by 
the much stronger field regions. The effect of this is that the 
largest scale field connects to the star at higher latitudes 
than would be expected when using a dipole magnetic field 
model. Thus by considering complex magnetic field topolo- 
gies we find that the fractional open flux through the stel- 
lar surface is less than would be expected from dipole field 
models. The reduction in the fractional open flux is mainly 
attributable to the larger closed flux for the complex fields, 



which arises due to the numerous small magnetic loops close 
to the stellar surface. For both V2129 Oph and BP Tau there 
is less total open flux than for a dipole field, however, the 
difference is small enough (see Fig. [5} that there is unlikely 
to be any significant implications for the amount of angular 
momentum tha t can be carried away from suc h systems by 
a stellar wind (JMatt fc Pudrit j l2005bl . l2008al lbh. However, 
more spectropolarimetric data for many different T Tauri 
stars will be required to confirm this. 

If the extent of stellar magnetospheres are limited to 
well within corotation then accretion may proceed along the 
more complex regions of the stellar surface field. The result 
of this would be the formation of many l ow-latitude accre- 
tion hotspots (|Gregorv et al.ll2005l . l2006aT ) . There is at least 
one star which shows evid ence for such equatorial spots, 
GQ Lup (JBroeg et al.l 120071 ). Ho wever, for V2129 Oph and 
BP Tau the spectroscopic data of lDonati et all i|2007l . l2008ah 
clearly showed that the bulk of accreting material is car- 
ried into high latitude regions of the ste llar surface. I n order 
to e xplain such high lati tude hotspots, IDonati et al.l (|2007l ) 
and lJardine et al.l (120081 ) demonstrated that for V2129 Oph 
the source surface must be set to at least 6.0 7?*, while 
for BP Ta u the stellar magnet osphere must extend to at 
least 4 7?* l|Donati et al.ll2008al ). Selecting the location of 
the source surface sets the structure of the stellar magnetic 
field, although the inner disc may still be truncated closer 
to the stellar surface. Using the assumption that inner disc 
is located at the point where the torque due to viscous pro- 
cesses in the disc is comparable to the magnetic torque due 
to the stellar magnetosphere we estimated the disc trunca- 
tion radius using the extrapolated fields. For BP Tau we 
found that the inner disc would be truncated close to the 
radius predicted for a dipole field, whereas for V2129 Oph, 
the inner disc is likely to be closer to the star. However, 
with large uncertainties in the assumed mass accretion rate 
(which influences the torque due to viscous processes in the 
disc), and with the static field structures considered here, 
such calculations are only of qualitative value. Further 3D 
MHD si mulations, such as those already performed by, for 
example I Long et al.l (|2008l ). are required in order to model 
how magnetic fields with an observed degree of complexity 
can disrupt and influence the structure of planet forming 
discs, and will represent a major future development in this 
field. 

Although the radius of the source surface is essentially a 
free parameter of our model, the observed hotspot locations 
provide constraints. Once the source surface has been se- 
lected, however, this then also sets the structure of the X-ray 
emitting surface field, allowing predictions t o be made re- 
gardi ng the stellar X-ray emission properties (|jardine et al.l 
120081 ). Such predictions will allow direct testing of our field 
extrapolation model with future independently obtained X- 
ray observations. A further test of the model will be the 
simulation of accretion related emission line profiles, and 
are currently being undertaken. 

V2129 Oph, despite its young age, is massive enough to 
have developed a radiative core. In contrast to this BP Tau 
is likely to be completely convective. Although the field of 
BP Tau is more complex than a dipole, it is much sim- 
pler than that of V2129 Oph. The differences in the field 
structure of both stars likely reflect the different internal 
structure and the process by which their magnetic fields are 
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generated and maintained. For stars with radiative cores 
the magnetic field is likely generated in the shear layer be- 
tween the core and the outer convective envelope. In the 
absence of such a shear layer, it is difficult to explain how 
fully convective stars can generate and maintain large-scale 
almost axisymmetric kilo-Gauss fields, alth ough recent theo - 
retical models suggest that this is possible (|Browning|l2008h . 
With only three magnetic surfa ce maps available to date , 
V2129 Oph and two of BP Tau (JDonati et alJl200ll2008al ). 
it is not possible to determine whether this is a general 
feature of all T Tauri stars, however it is si milar to what 
it foun d for low mass main sequence stars (JDonati et al.l 
|2008bJ). More magnetic surface maps of T Tauri stars with 
varying stellar parameters are now required to test these 
ideas fully, and are an essential requirement to advance 
our understanding of stellar magnetism on the pre-main se- 
quence. 
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